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O Optical dispersion spectra at energies up to 30 eV play a vital role in understanding the chirality-dependent van der Waals 
London dispersion interactions of single wall carbon nanotubes (SWCNTs). We use one-electron theory based calculations to 
obtain the band structures and the frequency dependent dielectric response function from 0-30 eV for 64 SWCNTs differing in 
radius, electronic structure classification, and geometry. The resulting optical dispersion properties can be categorized over three 
i "! distinct energy intervals (M, 7t, and a, respectively representing 0-0.1, 0.1-5, and 5-30 eV regions) and over radii above or below 
the zone-folding limit of 0.7 nm. While K peaks vary systematically with radius for a given electronic structure type, o peaks are 
f-j independent of tube radius above the zone folding limit and depend entirely on SWCNT geometry. We also observe the so-called 
metal paradox, where a SWCNT has a metallic band structure and continuous density of states through the Fermi level but still 
behaves optically like a material with a large optical band gap between M and K regions. This paradox appears to be unique to 
armchair and large diameter zigzag nanotubes. Based on these calculated one-electron dielectric response functions we compute 
. and review Van der Waals - London dispersion spectra, full spectral Hamaker coefficients, and van der Waals - London dispersion 
interaction energies for all calculated frequency dependent dielectric response functions. Our results are categorized using a new 
£h optical dielectric function classification scheme that groups the nanotubes according to observable trends and notable features 
I (e.g. the metal paradox ) in the 0-30 eV part of the optical dispersion spectra. While the trends in these spectra begin to break 
' ^ down at the zone folding diameter limit, the trends in the related van der Waals - London dispersion spectra tend to remain stable 
^ all the way down to the smallest single wall carbon nanotubes in a given class. 

l— — 1 1 Introduction their electronic structure and characterizafionESKl The rea- 

son is that SWCNTs are a unique class of materials in which 
p. Since their discovery in 1993 by the Iijima and the Bethune small changes in their atomic structure and geometry char- 
, — I groupsDEl single wall carbon nanotubes (SWCNTs) have re- acteristics (quantified by the chirality vector (n,m) or equiv- 
O ceived enthusiastic attention. Their intriguing mechanical and alently by [n,m, type], when we want to stress the type of 
electrical properties make them ideal for a wide range of appli- the nanotube [} determines whether they behave as a dielec- 
cations, from metal oxide semiconductor field effect transis- trie, a metallic conductor or a semiconductor. They are also 
i tors (MOSFETs) to novel drug delivery systems^ and have very narrow (typically < lnm wide) and inexpensive to pro- 
inspired strong scientific and industrial interestEEEl. The area duce, making them ideally suited for many nano-electronic 
that has arguably received the most attention and focus is applications while widespread commercial and industrial use 

has yet to be achieved. While there has been some recent 
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mentsEEKH s i v ing the assembly problem, which has been 
addressed on various levels beginning shortly after the discov- 
ery of CNTs^, presents an important next step based on the 
detailed molecular structure of various SWCNTs that deter- 
mines their optical properties and through them their interac- 
tions: a subject of this review. 

Understanding the optical dispersion properties of SWCNT 
is crucial to their role in several long-range interactions, espe- 
cially the van der Waals - London dispersion (vdW-Ld) inter- 
action. Each SWCNT is uniquely defined by its chirality vec- 
tor (n,m), which denotes the circumferential wrapping direc- 
tion along a graphene sheet containing sp 2 hybridized carbon 
bondslSEEZl Thus each (n,m) combination denotes a different 
geometrical construction that is responsible for the wide vari- 
ety of electronic structure, phonon, and Raman behavior 28 30 . 
The same diversity that makes them appealing also for perfor- 
mance leads to selectivity and specificity in sorting and place- 
ment. This is seen in AC/DC experiments 3 ^, dielectrophore- 
sis^2l, and anion IEC-SHSl Each chirality vector gives SWC- 
NTs a unique set of properties (e.g., band gap, van Hove singu- 
larities (vHs), etc.), but their common composition and chem- 
ical bonding coupled with symmetry considerations makes it 
possible to categorize them by their particular featureg 27 l 30 i 

A comprehensive study of CNT interactions can provide the 
fundamental conceptual framework required for proper under- 
standing of experiments on their separation, placement and as- 
sembly. Among these interactions, the universal and longest 
range, the van der Waals-London dispersion (vdW-Ld) inter- 
action, merits most detailed scrutiny. First, because these in- 
teractions are long range and ever-present in nature, always 
contributing to and influencing the overall interaction between 
charge-neutral bodies 33 , second, previous resultsEKD show 
that the magnitude of this interaction for various CNTs can 
vary by 20% even in simple systems, making it a convenient 
candidate to control experimentally, and third, SWCNTs have 
well known trends in their electronic structures as a function 
of chirality 19 , which could then act as the source for chirality- 
specific vdW-Ld interactions. 

During the last few years it has become possible to compute 
quantitative direction-dependent optical dispersion properties, 
such as the frequency dependent complex dielectric function 
e(co) = e'(co) + ie"(&), for a diverse set of isolated SWC- 
NTs, for which spectroscopic techniques even as powerful as 
electron energy loss spectroscopy (EELS)™, which requires 
bulk sample of identical SWCNTs in a bundled cannot be 
applied directly. The use of one-electron theory based calcu- 
lations allows us to circumvent these shortcomings but still 
use available experimental results to compare and confirm the 
integrity of our spectral data^^ES. We were able to resolve 
the axial and radial direction of SWCNTs to obtain the optical 
dispersion properties out to photon energies of at least 30 eV. 
Also, we were not limited to the SWCNTs that are easily pro- 



duced by existing preparation method s ! ' \ which allowed us 
to study a larger set of data than would otherwise be possible. 

Early attempts to obtain this information via ab initio codes 
produced unreliable results above 20 eV and were limited to a 
single nanotube 41 . Some computational studies were able to 
investigate a large quantity of chiralities, but they were con- 
fined to dispersion properties below 6 eV because they relied 
on the tight binding approximation, which artificially distorts 
the bands'^. New ab initio results are emerging that have ac- 
curate dispersion properties along the entire interv al, bu t con- 
tinue to be confined to a select number ofchiralitieslSEl. Here 
we calculate and review the optical dispersion properties, rep- 
resented as the imaginary part of the dielectric response func- 
tion, e"(K>), up to 30 eV for 64 SWCNTs that have varying di- 
ameter, electronic structure classification, and geometry. We 
follow trends and variations to give insight for system design. 
Later, we use these dielectric response functions to analyze 
consequences of different spectral properties in vdW-Ld inter- 
action energies. 

Trends and deviations from these trends in SWCNT optical 
dispersion spectra are categorized into a new SWCNT classifi- 
cation scheme, based on a combination of electronic structure, 
geometry and curvature. The links between chirality, cuttings 
lines, band diagrams, density of states and optical dispersion 
properties are extended to the magnitude of vdW-Ld interac- 
tion as quantified by the Hamaker coefficients for chirality- 
dependent vdW-Ld interactions. We highlight the key remain- 
ing features linking the dielectric function £"(co), with the con- 
sequent vdW-Ld interaction energy. 

2 Electronic structure and optical dispersion 
spectra 

2.1 Method 

The method for calculation of the electronic structure and 
optical properties is the first principles orthogonalized linear 
combination of atomic orbital (OLCAO) method which is de- 
scribed in detail in Ref. 4 ^. The OLCAO method is applied to 
the calculation of the electronic structure and optical disper- 
sion spectra of 64 SWCNTs. It yields the imaginary part of 
the frequency -dependent dielectric response function, e"(co), 
calculated within the random phase approximation (RPA) of 
the one-electron theory 47 . The optical transitions from the va- 
lence band (VB) to the conduction band (CB) are calculated 
from the ab initio wave functions using the extended basis 
(EB). The ab initio nature of the wave function ensures the 
strict adherence of the selection rules via the momentum ma- 
trix elements. If an optical transition is symmetry forbidden, 
it automatically has a zero momentum matrix and no transi- 
tion will occur even if such a transition appears to be available 
via the density of states (DOS). A key example that we will 
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show later is when a SWCNTs DOS has a continuous series 
of states across the Fermi level and yet behaves like a semicon- 
ductor with no transitions below 1 eV. All the re levant details 
of the calculation have been described i n IHEIIE] 

Though not taken into account by the one-electron method, 
we are fu lly aware of the importance of the excitonic ef- 
fects™ 

in the low dimensional systems such as graphene 
and They tend to introduce additional peaks in 

the energy range below 5 eV or slightly shift the position 
of some other peaks. From a pure theoretical point of view, 
the effect of excitons has been intensively studied and can be 
approached from many different angles including two-photon 
excitation, excitons of delocalized higher levels etc.. An im- 
portant question in this respect is to what extent the exci- 
tonic peaks can modify the optical absorption obtained from 
the one-electron theory? In some cases the traditional one- 
electron calculation of optical properties gives dielectric con- 
stants or the refractive indices in good agreement with ex- 
periments on many insulating materials. However, though 
the excitonic effects in the case of SWCNT optical spectra 
are important and eminently measurabl e ^ 28 ! 29 ! it is neverthe- 
less not clear how the changes in optical spectra due to exci- 
tonic effects will affect the overall Hamaker coefficients and 
the trends in the vdW-Ld interaction between different classes 
of SWCNTs. Since the vdW-Ld interaction is a functional 
(or more precisely a discrete sum over the Matsubara fre- 
quencies) of the frequency-dependent dielectric response, it 
is reasonable to expect that such effects will be small and that 
the strength of the interaction will depend more on the global 
properties of the dielectric response delimited by the sum rule 
then on local specifics, (see later results in section 6.3). 

The above conjecture is supported by the recent work of 
Hobbie and co-workers on empirical evaluation of attractive 
vdW potential in SWCNTs^. Even though they have also 
recognized the presence of excitonic effect in CNTs, the au- 
thors used a much simpler empirical approach by fitting the 
collective peaks to a superposition of individual Lorentzians. 
They have concluded that for the semiconducting nanotubes, 
neglecting the three optical resonances reduces the Hamaker 
coefficient by roughly 5%, while neglecting the Drude term 
decreases the Hamaker coefficient by roughly 2%. For the 
metallic SWCNTs, neglect of either of these terms reduces 
the Hamaker coefficient by roughly 3%. In this view, exci- 
tonic effects certainly have a small but measurable effect and a 
larger question is now how they influence the plasmons. From 
the ab initio perspective, a first-principles approach that cor- 
rectly accounts for electron-hole interactions in quasi-lD sys- 
tems would thus provide a more accurate computational foun- 
dation, but only to the extent that it also accurately describes 
the two plasmons. On the other hand, practical applications to 
real problems using optical spectra obtained from one-electron 
calculation is an attractive alternative as indicated in another 



recent paper on adhesion in silicon wafers^. 

The fundamental theory of the calculation method within 
the framework of the DFT is sound. Although many-body in- 
teractions are not considered and there may be concerns about 
the underestimation of the band gap, they are less important 
in this study since our main focus is on the vdW-Ld spectra 
of a large number of CNTs up to high frequency and their 
systematic correlations to the structure. In contrast, in the 
tight binding approximation (TBA) method where no accu- 
rate wave functions are involved, the band gaps are typically 
adjusted in order to best describe and represent the states very 
near the band gap. States several eV beyond the band gap 
are ignored or poorly represented^. In addition, TBA mod- 
els would not be able to predict novel departures from the 
main SWCNT structural trends. There are several examples 
of this, such as 1) the [5,0,s] actually is a metallic conduc- 
tor instead of being semiconductor-like^, and 2) the [6,0,m] 
has a large degree of hybridization of graphenes sp 2 and % 
bonds into sp 3 bonds and thus a band diagram unlike any other 
armchair SWCNT 3 -^. Ab initio codes capture these major de- 
partures from the trends without requiring ipso facto arbitrary 
adjustments. On the other hand, inclusion of excitonic correc- 
tions based on many-body theory is not feasible in the OLCAO 
method which uses local atomic orbitals for basis expansion 
rather than plane waves. It is interesting to note that empiri- 
cal TBA can also be refined to investigate the excitonic effect 
in SWCNT as has been demonstrated by the work of Jiang et 
alM 

The OLCAO method is an all-electron method with an eco- 
nomical use of the basis expansion. The full secular equation 
is diagonalized to obtain all electron states including high en- 
ergy CB states. If an iterative procedure is used for the elec- 
tron states, only a small number of CB state wave functions 
are usually obtained. To obtain higher energy CB states in an 
iterative approach for systems with a large number of atoms in 
the unit cell could be computationally prohibitive. The dipole 
matrix elements of transitions between VB and CB are ac- 
curately calculated from the ab initio wave functions and ex- 
plicitly included in the optical dispersion spectra calculation, 
which automatically impose the selection rules for the transi- 
tion. 

The OLCAO method is versatile and can be applied to al- 
most any material whether it is a metal or an insulator, a liquid 
or an amorphous solid, an inorganic crystal or a biomolecule, 
having an open structure or a compact one, or whether it con- 
tains light elements (such as H or Li) or heavy atoms (such 
as rare earth elements). The method avoids the use of atomic 
radii in the calculation, which could be problematic for com- 
plex materials with different local bonding. Finally, the use of 
local atomic orbitals facilitates the interpretation of the physi- 
cal properties where K or o states play a crucial role in the opti- 
cal spectra of SWCNTs. The efficient evaluation of multicen- 
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ter integrals in analytic forms makes the method highly effi- 
cient and eventually applicable to large complex systems such 
as DNA molecules in a fluid medium. For example, applica- 
tion to the calculation of large anisotropic optical absorption 
in herapathite crystal ((C2oH24N202H2)4 x C2H402 x 3SO4 x 
2I3 x 6H2O, a complex multi-component organic crystal with 
988 atoms in the unit cell yield excellent agreement with the 
measured data^l Application of the OLCAO method to the 
spectroscopic properties of a large model of super-cooled wa- 
ter is equally successful^. 

2.2 Case Study: SWCNT: [7,5,s] vs. [8,2,m] 

Each SWCNT chirality (n,m) defines a unique geometrical 
construction that results in a distinct set of electronic and 
optical dispersion properties. Because of their similarity to 
graphene in composition and structure, there are some fea- 
tures common among all SWCNTs, however, there are also 
two fundamental differences: circumferential periodicity and 
curvature. Periodicity along the circumference allows only 
discrete lines of states out of the continuum available for 
graphene, while curvature placed upon a flat graphene sheet 
re-hybridizes the sp 2 lattice to include sp 3 like states^. 

For an infinite diameter SWCNT, the number of available 
states is a continuum equivalent to nearly flat sp 2 sheets of 
graphene. The only differentiating feature among various in- 
finite diameter tubes would be the bonding direction along 
the axial direction. For a large, but finite diameter SWCNT, 
the effects of periodicity would cause a set of discrete al- 
lowable states in the axial direction. However, because the 
carbon bonds would still be relatively flat (i.e., still nearly 
perfect sp 2 bonds), it is possible to derive a SWCNTs band 
structure by selecting the allowable states from the band dia- 
gram of graphene. Such an approximation is called the zone- 
folding approximation; it is valid only for SWCNTs of diam- 
eter greater than 1 nm. Below this threshold, the curvature of 
the SWCNT wall becomes so pronounced that the bonds gain 
a significant sp 3 characteristic. For very s mall tubes , using a 
zone folding procedure can be problematic? 1 9 I 3 ° I 45 I 7 ° I. 

Each SWCNT chirality represents a unique geometrical 
construction of a ID tube of sp 2 carbon bonds created from 
a rolled up sheet of graphene. What makes SWCNTs such 
a highly studied set of materials is the number of material 
properties that significantly change as a function of chirality. 
Perhaps the best well known is metallic versus semiconduct- 
ing behavior, depending whether or not (n — m)/3 is an in- 
teger. Variations in properties propagate through several fea- 
tures, for example: linkages between chirality, cuttings lines, 
band diagrams, density of states, optical dispersion proper- 
ties, Hamaker coefficients, extending to chirality-dependent 
vdW-Ld interactions. Figure [T] shows these differences for the 
metallic [8,2,m] and the semiconducting [7,5, s]. 



While some of these relationships (e.g., chirality to elec- 
tronic band diagram) have been extensively studied and some 
predictive rules are widely used, much information at the op- 
tical level cannot be easily predicted a priori or explained sim- 
ply by knowing (n,m). Perhaps the most notable example is 
between the electronic structure and optical levels, which ex- 
hibit many non-linear relationships and other features that of- 
fer somewhat surprising results. The need is to determine the 
optical dispersion properties for most chiralities rather than 
approximate/derive them from a SWCNT of seemingly simi- 
lar characteristics. 

The introduction of ab initio e"(co) optical dispersion prop- 
erties enables such an analysis, which would be very difficult 
experimentally. Fortunately, the experimental data that do ex- 
ist (notably the EELS data by Stephan et al.^j) compare favor- 
ably to our results. Although they specify tubes by diameter 
and not chirality, they do note major peaks that appear consis- 
tently across tubes of all diameters. Most notably, their results 
show three major peak areas: band gap, transitions at 4 eV, 
and transitions at 15 eV. There are also low flat areas between 
5-12 eV and beyond 20 eV. These results are in qualitative 
agreement with our [7,5, s] and [8,2,m] ab initio optical dis- 
persion properties. The only major difference appears to be in 
the balance of the e"(co) peak heights and total areas within 
the 7t and o regions. However, these differences may resolve 
when experimental data can be more readily divided into the 
radial and axial components, or perhaps a finer experimental 
energy resolution would help. 

2.3 SWCNT classification schemes 

Table 1 lists the 64 SWCNTs studied here, with their fig- 
ures of merit (n, m, radius, angle, geometry, and number of 
atoms per lattice repeat), and the three classification schemes 
(zone folding, Lambin 19 , and Saito 9 ). Our proposed disper- 
sion classification (not in Table 1) is essentially a combina- 
tion of the Lambin and geometrical descriptors. Although this 
technically leads to nine possible combinations, only five of 
them are allowed for smaller diameter SWCNTs. For instance, 
an armchair SWCNTs electronic structure is never semicon- 
ducting unless there are added defects, functional groups, 
or other influential external source such as a strong external 
field™!. For the larger-diameter limit, the divisions may 
become redundant as the semimetal gap becomes so small 
(i.e., less than 0.01 eV) that it essentially blurs the line of the 
metal/semimetal division. 

Figure [T] shows the cutting lines for the SWCNTs repre- 
senting an example of two types of the dispersion classifica- 
tion, [7,5, s] and [8, 2, m]. It is the relationship of cutting 
line angle, location, and density that determines if and when 
the cutting lines cross/graze the K points at the edge of the 
BZ and result in a metal/semimetal electronic structure. Nan- 
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Energy (eVJ Energy (eV) 

Fig. 1 (color on line) Layers of abstraction of SWCNT properties ranging from [n,m] to the optical dispersion spectrum e"(co) for the [7,5,s] 
semiconductor and [8,2,m] metal, (a) The geometrical construction, (b) A chirality map showing all 5 types of SWCNT. (c) The cutting lines 
across the Brillouin zones for two different types of SWCNT, [7,5, s] and [8, 2, m], within the dispersion classification. The dot represents the 
K point in the Brillouin zone, (d) The optical dispersion spectrum e"(co) in both the radial and axial direction. Frequency in eV. 



otubes of very different cutting line density and angle (i.e., 
the [6,5, s] and the [9,l,s]) can have different cutting line an- 
gles and yet similar band gaps. Conversely, tubes that appear 
to have very similar cutting line angles and density (i.e., the 
[9,0,m] versus the [10,0,s]) fall into different classifications 
altogether. Such is the delicate nature of the electronic struc- 
tures dependence on (n,m). The cutting line representation in 
the BZ has served as the basis for every major SWCNT classi- 
fication thus far (i.e., simple, Lambin, and Saito). In the zone- 
folding classification, the effects of curvature on the carbon 
sp 2 bonds are ignored and the only thing that differentiates the 
basic metallic and semiconducting classifications is whether a 
K point is crossed. This results in the simple (n-m)/3 = integer 
rule^l The Lambin system acknowledges the changes in sp 2 
bond angles and lengths as a result of the cylindrical nature 
of SWCNTfPH. This curvature causes a slight distortion in 



the structure and the resulting BZ contains cutting lines that 
barely miss or graze the K-points and dont directly crossing 
them. This opens up the tiniest of band gaps (ranging from 
0.02 to 0.05 eV) for tubes that the zone-folding classification 
would predict as perfect metals with no band gaps. Figure 
[2] illustrates the cutting line diagrams and resulting electronic 
band structures of these three types. The [9,3,m] has a barely 
visible 0.02 eV band gap. For comparison, the Lambin model 
predicts a gap of approximately 0.07 eV. 

Finally, the Saito method 9 returns to the zone-folding 
scheme (i.e., ignores the effects of the structural relaxation), 
but it considers additional symmetry effects that can further 
differentiate tubes into numerous grouping For example, 
the semiconductors can be subdivided into two classes based 
on the relationships of the first cutting lines position relative to 
the K point. Neither the proposed dispersion classification nor 



This journal is ©The Royal Society of Chemistry [year] 



Journal Name, 2010, [vol], 1-J22] | 5 




classification system. Note the subtle difference between the [6,6,m] 
and the [9,3,m]. The [6,6,m] has no band gap whereas the [9,3,m] 
has CB and VB just miss touching. Thus the [9,3,m] is actually a 
small band gap semiconductor . 



the available classifications can adequately predict the disper- 
sion properties of the smallest SWCNTs. At this limit, there is 
a strong reintroduction of sp 3 -like bonding, which can engen- 
der new bands within the electronic structure as well as signif- 
icantly shift other features. For this reason, [5,0,s] is metallic 
despite no classification predicting it 45 . It appears, however, 
that there are very few of these exceptions. Therefore, one 
should always seek additional experimental or ab initio con- 
firmation when describing the dispersion properties of very 
small SWCNTs. 

2.4 Influence of structural parameters on optical transi- 
tions 

Despite the shortcomings in each of the given classifications, 
they are powerful tools that can help to map out and under- 
stand global trends. The first trend is that the packing density 
of cutting lines varies systematically as a function of radius, 
creating a diameter-dependent band gap and trends in the ax- 
ial e"(co) peaks. Figure [3] shows e"(co) for three different arm- 
chair SWCNTs, [4,4,m], [14,14,m] and [24,24,m]. Increas- 
ing the SWCNT radius allows more lattice translations in the 
circumferential direction and thus a larger degree of cutting 
line packing. Dresselhaus and Saito quantified the density of 
cutting 928 and showed that this increase in the cutting line 
density means that the cutting line closest to the K point will 
push closer to it. Therefore, the larger the SWCNT radius, the 
larger its cutting line density and the smaller the difference 
between its first conduction and valence bands. This causes a 
systematic shift in th e DOS and e"(co) peaks to lower energy 
(a well known effect)^ 2 ^^. Additionally, it means that larger 
tubes will have more of these peaks because of an increased 
availability of cutting line transitions. Note that these shifts of 
e"(co) peak position as a function of diameter are only valid 
for the axial direction K peaks. In Figure[3ja), one can see that 
aU the SWCNTs have a final e"(co) peak (the n peak) at 4.20 
eV in the 71 transition region. This is due to the cutting lines 
reaching some common border for all armchair SWCNTs. All 
peaks beyond this 4.20 eV Kq peak, in the 5-30 eV a transition 
region in Figure|3](b) are invariant with change in radius. 

Before investigating the differences in these two regimes, 
first consider the known variation in £"(co) peak locations be- 
tween the metal and semiconductor classifications in the K 
regime. When a K point is crossed, the next nearest cutting 
line is a full spacing away. Conversely, the semiconducting 
tubes have their first two cutting lines at 1/3 and 2/3 of this 
distance. The implications of this are: Metals cross a K-point 
where the CB and VB meet, allowing for the possibility of a 
finite-height metal peak all the way down to 0.00 eV. Sym- 
metry effects are the only thing that prevent it from appear- 
ing in the e"(co) optical dispersion properties. Additionally, 
metals have their next nearest cutting line at an energy three 
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(a) (b) 

Fig. 3 (color on line) The axial component of the optical dispersion spectrum e"(co). (a) e"(co) in the 0-5 eV region, (b) e"(oo) in the 0-30 eV 
region, both for the case of [18,18,m], [21,21,m] and [24,24,m] SWCNTs. Note that as the cutting line density increases for larger tubes, there 
is a systematic shift of the first £"(&>) 7Ii peak towards OeV. However, the optical dispersion spectrum £"(u)) above 5 eV appear to be largely 
invariant with respect to changes in diameter for large tubes. Frequency in eV. 



times larger than that of a semiconductor of equivalent radius. 
Therefore, while the K peaks in e"(co) from 0-4 eV will shift 
systematically as a function of radius for tubes of similar elec- 
tronic structure (i.e., metals or semiconductors), there is a big 
shift in the peak locations between the classifications'^^. 

As seen in Figure|3jb), the o transitions seen in the 5-30 eV 
optical dispersion properties seem to be invariant with respect 
to changes of CNT radius. But are the data in this regime the 
same for all SWCNTs or is there another parameter that can 
make a difference? If we compare a zigzag semimetal with a 
zigzag semiconductor, we have essentially fixed the geometry 
and tested whether there is a correlation with the electronic 
structure properties of the K regime. Likewise, if we compare 
zigzag semimetals with armchair metals, we have essentially 
fixed the electronic structure while varying the geometry H 
Comparing all three would allow us to check if there is a si- 
multaneous correlation with both geometry and 0-5 eV elec- 
tronic structure. Figure [3] (a) and (b) combined with Figure [4] 
(a) and (b) show just such a comparison of the three SWCNT 
types for both the 0-5 and 5-30 eV ranges. The % regime be- 
haves as expected with its dependencies on radius and zone- 
folding classification. The e"(co) peaks of the metallic and 
semiconducting zigzag SWCNTs are visually identical in the 
a regime, eliminating any correlation to the zone-folding clas- 
sification. However, there is a clear difference in e"(co) peak 
magnitude and position between the armchair and zigzag ge- 
ometries, suggesting geometry as the influential parameter in 



t Technically, semimetals and metals have differences in their properties at the 
Fermi level. However, we are more interested in whether the relationships of 
the cutting lines with the K-points also impact the 5-30 eV properties in any 
significant way, making the crossing/grazing differentiation irrelevant. 



the o regime. 

3 Features and analysis of spectral properties 

Because we have observed trends in spectral properties for the 
largest diameter SWCNTs (i.e., in the zone folding regime), it 
is instructive to ask how resilient they are. Do they work down 
to the smallest radius SWCNTs? Is there a particular cutoff 
for the range/diameter? Are the departures from the trends 
subtle or are they large, leading to a seemingly random set 
of properties among the chiralities? Figure [5] (a) to (d) shows 
e" (co) peak position vs. radius for the armchair SWCNTs from 
[3,3,m] to [24,24,m] in both axial and radial directions in two 
energy ranges separated at the critical radius of 0.7 nm. In the 
0-5 eV range there is a smooth trend down to the tiniest of 
SWCNTs. However, in the a regime, we see a large distortion 
of peaks for radii of less than 0.7 nm. New peaks appear to 
form and others appear to shift randomly^] 

In the case of SWCNTs we detect additional K — a* and 
a K* -like transitions due to curvature causing these other- 
wise independent directions to overlap. As a result, between 
the limit of 3.8 eV (all 71.-71* transitions below) and 10.0 eV 
(significant a — a* transitions above), it is not clear which 
particular transitions contribute most to the optical dispersion 
properties. The goal here, however, is simply to cathegorize 
the two distinctly different regimes found in Figs. [4] and|5] Ul- 
timately our selection of the 5 eV as the n/a cutoff was based 
on three major factors: 1) Large diameter SWCNTs typically 



% Note that a change in quantity of peaks doesnt necessarily change the overall 
magnitude. It simply is a useful parameter to monitor when studying changes 
among the various chiralities. 
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Fig. 4 Comparison of e"(co) optical dispersion spectrum in the (a) 71 ( 0-5 eV) and (b) a (5-30 eV) regions along the axial direction for 
metallic and semiconducting zigzag SWCNTs. Frequency in eV. 



have no major £"(co) contributions between 5-10 eV for the 
axial and/or radial directions. Thus we could arbitrarily se- 
lect any location along this interval without changing around 
the distributions of peaks contained in each of the groups. 2) 
Since the edge between K — K* and non K — K* related e"(co) 
contributions is the only hard edge known with certainty, it 
makes sense to select the lowest energy in the 5-10 eV in- 
terval. Therefore, our n classification for SWCNTs is purely 
71 — 71* based while anything above 5 eV is likely to have con- 
tributions from all 4 possibilities (71 — 71*, a — 71*, K — a*, and 
O - a*). 3) In Figure|5] the smallest SWCNTs have the largest 
curvature and therefore the highest degree of K — a* and a — K* 
overlap. This is likely the reason for the increase in e"(co) 
peaks appearing in the 5-10 eV limit. 

3.1 Metal Peak Analysis 

An ideal metal is characterized by a simple Drude model cor- 
responding to 

„ 4ti Ne 2 x 2 x 

e" (<*>)= j-=(0^— ^— j-. (1) 

m co(l+x 2 co 2 ) p co(l+x 2 co 2 ) 

where N is the electron density, m the electron mass, and co p 
the characteristic plasma frequency. This expresses the diver- 
gence of the dielectric response at vanishing frequencies. This 
feature, while not critical to our studies of SWCNTs, is a topic 
of interest in the formulation of electrodynamic forces due to 
the importance of conductivity affectsL^H^H i n computation, 
we can only approach the static limit asymptotically. 

In computation we can distinguish three types of behavior 
in this metallic peak energy range. We refer to these behaviors 
as null, Mo, and Mi. Null represents a SWCNT with at least 
a 0.1 eV optical band gap. For a Mo peak, e"(co) rises asymp- 
totically as the energy approaches eV, while a Mi peak will 



rise in a similar manner but then falls back to a value of e"(co) 
= at an energy in the range from 0.02 to 0.05 eV. It should 
be noted that very small band gaps (less than 0.05 eV) do exist 
in materials and in this work the Mi type metal peaks are not 
likely due to a computational artifact since we calculate the 
optical transitions and dipole matrix elements down to 0.01 
eV transition energy. 

Consider the optical dispersion properties of three SWC- 
NTs representing the armchair metal, zigzag metal, and chiral 
metal classes. In a zone-folding classification, these would all 
be considered metals and therefore the eV behavior should 
be indistinguishable for all three cases. The Lambin classifi- 
cation would state that the [9,9,m] armchair is a true metal, 
but the [18,0,m] and [9,3,m] would experience a slight gap of 
less than 0. 1 eV. Therefore we would expect to potentially see 
two types of behavior at eV. In actuality, we see all three 
metal behaviors: null (armchair), Mi (chiral metal), and Mo 
(zigzag metal). The most surprising observation is that of the 
[9,9,m] armchair. Despite having no band gap and a tiny but 
continuous DOS around the Fermi level, it exhibits no optical 
transitions from to 1.5 eV. This effect is somewhat paradox- 
ical because the armchair SWCNTs represent the only class 
with no band gap (not even a very small band gap like the 
zigzag and chiral semi-metals), and yet they exhibit optical 
band gaps as large as 1 eV. Also unexpected is the Mo peak 
of the [18,0,m] because zigzag metals are supposed to have 
a very small, but finite band gap and thus we would expect 
an Mi type peak. Perhaps structural relaxation would resolve 
this inconsistency and uncover an Mi type peak instead^. If, 
however, this proves to be a real and accurate result, then di- 
viding the metals into these three metal categories is justified 
because it will result in three types of vdW-Ld interactions. 

Because physical properties depend strongly on the exis- 
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Fig. 5 Plots of peaks in optical dispersion spectrum £"(&>) vs. the radius of SWCNTs. (a) axial direction from 0-5 eV in the 7t region. Smooth 
trends observed for all radii; (b) radial direction from 10-20 eV in the a region. No trends for small radius tubes; (c) axial direction from 0-5 
eV in the n region. Smooth trends observed for all radii; and (d) radial direction from 10-20 eV. No trends for small radius tubes. Frequency in 
eV. 



tence and magnitude of an Mo or Mi peak, we wanted to 
determine if such peaks scale with SWCNT radius. Figure 
|6^a) compares the e"(co) dispersion properties in the 0-1.5 eV 
range across 8 different zigzag metal SWCNTs. There is a 
clearly inverse relationship between the metal peak magni- 
tude and radius with the [9,0,m] having a Mo peak of approxi- 
mately 20,000. This magnitude systematically diminishes un- 
til [24,0,m], where the Mo peak opens up to a very weak Mi 
peak. For larger tubes, this peak disappears altogether. 

Do chiral metals experience a similar Mo/Mi trend as a 
function of radius? Not according to the metal peaks for 
the nine chiral metal SWCNTs listed in Table 1. The largest 
M peaks come from the [5,2,m] and the [10,4,m] SWCNTs, 
which respectively have the smallest and second largest radius 
within the group. The other seven chiral metals have Mi and 
Mo peaks that are weaker by 1 to 2 orders of magnitude. More 
large-diameter chiral metals are need to be studied, but they 
are more difficult to study as they often require up to an order 
of magnitude more atoms within the unit cell of the ab initio 
calculation. 



A Drude metal shows a 1 /co divergence of the dielectric re- 
sponse function (Eq. [TJ, giving a slope of -1 in a log-log plot 
of e"(fi)). Figure |6jb) shows this plot for five zigzag metal- 
lic single- wall CNTs. Instead of observing a slope of -1 (not 
shown), we observe a slope of -2. The significance of this 
difference is not clear . 

3.2 7i and o region analysis 

All large diameter SWCNTs show a TCo peak in e"(a>) at 4.0- 

4.3 eV that is invariant with respect to changes in radius for 
a given SWCNT class and direction (i.e. radial or axial). For 
armchair SWCNTs [see Figure |3ja) and (b)] the 7to peaks oc- 
cur at 4.19 eV and 4.21 eV for the axial and radial directions, 
respectively. However, for armchair tubes with diameters be- 
low the zone-folding limit, these peaks shift noticeably with 
the axial (radial) direction shifts slightly lower (higher) as a 
function of 1/radius. Why these peaks are pinned at the large 
diameter limit and shift in opposite directions at the small di- 
ameter limit is unclear. It is conceivable that all SWCNTs ex- 
perience some common edge in the BZ that begins to dissolve 
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(a) (b) 

Fig. 6 (color on line) (a) e"(oo) optical dispersion spectrum around eV for eight zigzag SWCNTs spanning a large range of diameters, (b) A 
log-log plot of e"(oo) versus energy for five zigzag metals. The data below 0.1 eV begins to express a flat slope of -2 , indicating a divergent 
behavior down to eV. Frequency in eV. 



as the cutting line density decreases and the bond distortion 
increases. 

The next important peak is associated with the optical band 
gap, labeled tc i . All other peaks between K\ and Kq are incre- 
mentally labeled K2, 713, etc., if they exist, [see Figure[3jb)]. 
The largest SWCNTs in this study have up to six K peaks in 
the axial direction. The radial direction typically has only Kq 
and fti peaks with relatively flat, low regions in between. 

Zigzag metals and semiconductors have a little more vari- 
ation in their Kq peaks as shown in Figure |4ja) and Figure 
|4|b). The values typically oscillate between 4.00 and 4.30 
for the largest zigzag SWCNTs in the axial direction, with 
the variation increasing versus 1/radius. Just like the arm- 
chair SWCNTs, the peaks in the radial and the axial directions 
shift in the opposite direction as the radius decreases. The 
Kq peaks for chiral metals and semiconductors are a bit var- 
ied. We have no data for SWCNT with radius greater than 0.7 
nm, or larger than the zone folding limit. It is unclear if they 
will have a stable 7to peak for a certain diameter range. As 
a whole, tc 1 — 715 peaks trend as expected, particularly above 
the zone-folding limit. Below this limit, some e"(co) peaks 
(most notably for the chiral geometries) distort and shift in 
unpredictable ways, underscoring the need to obtain this in- 
formation for each SWCNT rather than approximate it using 
the data of another chirality. In contrast, the armchair class 
of SWCNTs has 71 1 — K$ peaks that vary smoothly down to the 
smallest diameters. The similar geometry, electronic structure, 
and symmetry probably contribute to this resiliency. 

Like TCo, the a peaks are typically invariant as a function 
of radius within a given SWCNT class for diameters above 
the zone-folding threshold. Up to three o peaks are typically 
identified. The Gq peak tends to be around 14.5 with a sister 
a 1 peak around 13 eV. Armchair SWCNT geometries have an- 



other more rounded peak around 27 eV whereas the zigzag ge- 
ometries do not. The remaining peaks that arise for very small 
diameters are not readily named with this system because they 
can be sharp, often contain only a limited area under the curve, 
and are difficult to compare tube to tube. With a coarser res- 
olution in the e"(a>) data (e.g., 0.05 eV instead of 0.01 eV), 
many of these peaks simply disappear. We collectively name 
these a* because their significance depends on the context of 
the question to be answered. 

For the o peaks, the armchair SWCNTs have values of Gq 
and ai of 14.56 and 13.45 (14.56 and 12.69) respectively in 
the axial (radial) direction. Thus the location of the Oq peak 
is the same for both directions and the Gq shifts roughly 0.8 
eV. For the zigzag chirality (metal and semiconducting), the 
Gq peaks location is also at 14.56 for both directions. But 
the Gi peaks correspondingly shift about 0.1 eV to 13.60 and 
12.83 for the axial and radial directions, respectively. These 
shifts may seem small, but changes in position for peaks hav- 
ing large areas under the curve could have a large impact on 
vdW-Ld interaction. 

Thus far, our primary focus has been on trends of the dis- 
persion properties of a large population of categorized SWC- 
NTs. At the small diameter limit, there are no stable jiq or 
Gq peaks, and the rest of the peaks can vary in unpredictable 
ways. Increased curvature introduces more sp 3 hybridization, 
causing the electronic structure to depart strongly from that of 
grapheneESl. Also, the large decrease in cutting line density 
means that the allowed states within a given BZ are quite dif- 
ferent. Large optical dispersion variability is found when the 
different electronic structures are combined with the different 
allowable samplings. This is illustrated in Figure [7] which 
shows the e"(co) spectra of four very small diameter chiral 
tubes in the K and o regions. 
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(a) (b) 

Fig. 7 (color on line) e"(co) optical dispersion spectrum in the (a) 0-5 and (b) 10-20 eV range for four very small, chiral SWCNTs. Frequency 
in eV. 



4 General observations on the optical proper- 
ties of SWCNT 

Perhaps the most important contribution from the above anal- 
ysis is the quantification of the variations of e"(co) in the 0.0- 
0.1, 0.1-5, and 5-30 eV ranges. The basic trends that occur 
in the 0-5 eV range are established in the literature. The op- 
tical dispersion properties in this range are mostly correlated 
with electronic structure, which depends on the cutting line 
position (i.e., whether the K-points are crossed) and the pack- 
ing density (which shifts the peaks toward or away from eV 
based on the distance of the cutting lines from the K-points). 
The only major departure from the previous electronic stud- 
ies is that of the Mo regime. Here we discovered that some 
metals (notably armchair SWCNTs) do not exhibit a standard 
metal peak despite having the available valence and conduc- 
tion band states. This results in the metal paradox where the 
armchair metals are optically like semiconductors^! 

In contrast to the M and K regimes, the optical properties in 
the 5-30 eV range are largely dictated by the effect of under- 
lying geometry on band to band transitions at higher energy. 
This means all zigzag SWCNTs at the large diameter limit 
will have identical properties above 5 eV despite having dis- 
tinct behaviors (i.e., null, Mi, or Mo metal peaks and metal 
versus semiconducting K peaks) in the lower energy regimes. 

Having extended our focus beyond DOS and electronic 
structure, the interesting symmetry considerations outlined by 
Dresselhaus, Saito et al. do not necessarily manifest them- 
selves as notable differences in the optical properties 9 . Certain 
optical transitions are symmetry forbidden, unlike the analy- 
sis of Katura plots and the VB to CB transition tracking with 
radius^. This is a hazard of stopping at the DOS level anal- 
ysis and approximating the optical properties without explic- 



itly calculating the dipole matrix elements. For example, the 
known differences (metal versus semiconducting position of 
7ti) and trends (7ti — 715 shifting with radius) lead to easily 
observable differences in key optical properties, such as the 
index of refraction. In the small diameter limit (where the 
aforementioned trends have broken down) SWCNTs like the 
[6,2,s], [6,l,s], and [5,l,s] have little correlation in the e"(co) 
peaks leading to a distinct set of properties for each of these 
smaller tubes. Such variations are important for material se- 
lection and experimental design. While we cannot make broad 
predictions without performing calculations, the calculations 
for vdW-Ld interactions should be straightforward and can be 
used to find the relevant variations (see below). 

Our new classification is different from the existing ones, 
both in terms of the particular categories selected and the rea- 
sons for their selection. For example: Lambin et al., correctly 
moved beyond the simple zone-folding classification by iden- 
tifying the armchair metals as true metals while ascribing very 
small gap semiconductor characteristics to chiral and zigzag 
metals 19 . But it is the armchair that behaves optically like 
1+eV semiconductors. The zigzag and chiral semimetals with 
their sub 0.05 eV band gaps exhibit some form of a metal 
peak. Thus the division between armchair and non-armchair 
SWCNTs is correct, but the correlation is opposite from ex- 
pectation. The Dresselhaus/Saito system, despite having in- 
teresting geometry and symmetry considerations relating di- 
rectly to (n,m), appears to be less accurate as a categorization 
tool for the optical properties in the range. As an example, 
the M-2p metal category in this classification contains both 
armchairs and chiral metals. But as seen earlier, the metal 
peaks of these two groups are different and will lead to differ- 
ent vdW-Ld interactions. This discrepancy is a result of the 
zone-folding approximation. Despite these shortcomings in 
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explaining the e"(co) optical dispersion properties, the descrip- 
tors used to create this classification are useful in describing 
more quantitative differences in cutting line position, angle, 
and packing density. Whether some e"(co) optical dispersion 
properties link with the symmetry differentiation contained 
with this classification can only be determined when results 
on more chiral semimetals and semiconductors are available. 
Until then, geometry is still a far more important parameter 
for the 0-30 eV optical properties^. 

5 van der Waals - London dispersion interac- 
tions 

Represented by the chirality vector (n,m), chiral properties of 
SWCNTs are connected with their vdW-Ld interactions that 
depend on their optical properties through the dielectric func- 
tion, e"(co). The Rramers-Kronig transform of e"(co) leads 
to the vdW-Ld spectrum (vdW- Lds) e(it,), i.e., the dielectric 
response function at the imaginary values of the frequency 
that serves as the basic input for the Lifshitz theory of vdW- 
Ld interactions! 34 ! 34 ! ^^. Above, we addressed and discussed 
in detail the connection between the chirality vector (n,m) and 
the optical dispersion spectra e"(co). While non-intuitive fea- 
tures can occur when one tries to categorize the S WCNT prop- 
erties based on their electronic spectra (e.g. the metal paradox 
for armchair SWCNTs), most of the relationships between the 
S WCNT atomic structure and geometry and their optical prop- 
erties are conceptually straightforward. We now proceed from 
electronic structure and optical properties of SWCNTs to their 
vdW- Ld spectra and the ensuing vdW- Ld interactions con- 
necting chiral properties of SWCNTs with long-range interac- 
tions. 

5.1 van der Waals - London dispersion spectra 

Connecting the imaginary part of the dielectric function e"(co) 
with the corresponding vdW-Ld spectrum, e(?£), the Rramers- 
Kronig (RR) relation^ 

7t Jo co- + q z 

describes the connection between material's response to exter- 
nal electromagnetic fields of a given frequency and the vdW- 
Ld spectrum, which characterizes the magnitude of sponta- 
neous electromagnetic fluctuations at the frequency t,. In gen- 
eral e(it,) is a real, monotonically decaying function of the 
imaginary frequency i^S. 

Note that the integration in Eq. [2] is over an infinite fre- 
quency range. In practice, it is usually impossible to know the 
frequency response of a material accurately over the whole 
frequency range. Fortunately this is unnecessary as long as all 



the inter-band transition energies determining the optical dis- 
persion spectrum are either known or properly approximated. 
The e"(co) at high energies in the optical domain needs to be 
pronounced in order to contribute to the resulting vdW-Lds in 
any significant way because of the heavy damping in the de- 
nominator of Eq. [2] The SWCNT e"(oa) spectra above 30 eV 
described in this paper were relatively small in magnitude and 
noisy. Therefore a cutoff energy of 30 eV was selected to best 
achieve the desired accuracy when evaluating the Rramers- 
Rronig transforms. 

There exist distinct conterbalancing effects between high 
energy wide peaks or low energy narrow peaks, and because 
e"(co) and b(i£) are connected via an integral transform, it is 
possible to have identical vdW-Ld spectra with very different 
dielectric functions. Furthermore, excitonic effects that mod- 
ify e"(co) with peak structures 50 will have in general only a 
small effect on e(;^) 66 , which remains a smooth monotoni- 
cally decreasing function of its argument. In fact, as is clear 
from Eq. [2] every peak will contribute to the vdW-Ld trans- 
form approximately additively, as the area under the peak, 
scaling inversely with i.e. inversely with the square of 
the Matsubara frequency. Since it is the Matsubara frequency 
summation in the optical regime of the vdW-Ld transform, that 
determines the Hamaker coefficients and the strength of vdW 
interactions, see Eq. [3] below, the corresponding excitonic 
effects on these, while of course present, will be by neces- 
sity small. This statement is substantiated by a more detailed 
quantitative analysis based on full measured SWCNT optical 
spectra containing an explicit excitonic partial 

5.2 Hamaker coefficient 

The Hamaker coefficient measures the strength of the vdW-Ld 
interaction energy between two interacting materials'^EMI] 
It is defined with respect to the geometry or configuration of 
the interaction setup and the shape of the interacting bodies. 
Its magnitude depends upon the optical contrast between the 
interacting materials and the bathing medium. It varies as a 
function of separation, a consequence of retardation. In the 
limit of small separations it can be spoken of as a a constant. In 
what follows we shall consider only the vdW - Ld interactions 
at such small separations, either between two parallel SWC- 
NTs or between a single SWCNT parallel to an anisotropic 
planar substrate 3 -40fil 

The cases of larger and very large separations between two 
SWCNTs and between a SWCNT and a planar substrate have 
been dealt with in our other public ations that also contain full 
derivations for all separation limits 34 l 34 l 36 l 82 l. Our analysis is 
strictly applicable only in the limit of infinitely long SWC- 
NTs; for a consistent analysis of finite size effects and ideal 
metallic static dielectric response one needs to consider an ex- 
act general multiple scattering formulation of the vdW - Ld 
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interactions, such as derived by Emig and coworkers^. 

Consider first the vdW - Ld interaction free energies be- 
tween two identical infinitely long parallel SWCNTs at small 
separations where retardation effects can be neglected^. We 
define the dispersion anisotropy parameter for a cylindrical 
SWCNT, f (co), as 



e c n(^)-e c ± (^) 



(3) 



where e c || ) j_(j^) are the parallel (with respect to the cylinder 
axis) and the perpendicular components of the vdW - Ld spec- 
trum of the cylinder. The van der Waals interaction free energy 
per unit length between two such parallel SWCNT of radius 
a, at a separation £ <C a, has the fomJS] 



g(l,B;a) = - 



24 p/ 2 



(4) 



where the sum of coefficients J?(°) = \k B T rf® and ^ = 
\k B T is 

2k ^ Jo 



with 



£' ± (4,) V 1 + f (l%„ ) COS 2 \|/ + 

Here e„,(i^„) is the vdW - Ld spectrum of the intervening 
medium. The vdW-Ld response functions in above equa- 
tions are given at a discrete set of Matsubara frequencies 
\n = B % " » where k B isthe Boltzmann constant and h is the 
Planck constant divided by 271. At room temperature, the in- 
terval between neighboring Matsubara frequencies is approxi- 
mately 0.16 eV. The prime in the summation signifies that the 
first, n = 0, term is taken with weight I . 

In general at a finite mutual angle of the dispersion 
anisotropy axes 8 the vdW interactions are codified by two 
values of the Hamaker coefficient, i.e. J4.^ and M 2 \ whereas 
in the par allel geom etry they are given by a single value 

Sl <0) I J^ 2 H34I34I36I82I 

The most important aspect of the connection between the 
Hamaker coefficients and the magnitude of the vdW-Ld in- 
teractions is the optical contrast, or dispersion mismatch 
A (e c ± , £„, , Y ') in Eq. |6] essentialy proportional to the differ- 
ence of the dielectric response function over their sum at ev- 
ery abrupt dielectric interface. The Hamaker coefficients for 
anisotropic materials are thus proportional to the angular \\t 
integral of the product of all the dispersion mismatches in the 
system. Clearly, if two materials at each dielectric interface at 
zero mutual angle have identical dielectric responses at a given 



Matsubara frequency, the mismatch term goes to and there 
is no contribution to the Hamaker coefficient at that particular 
frequency. Thus the magnitude of the Hamaker coefficient is 
dictated by the degree of dispersion contrast at both interfaces. 
On the other hand, the sign of the interaction (i.e. attraction 
vs. repulsion) is determined by the sign of the dispersion mis- 
match at the two interfaces if they are not identical. Depend- 
ing on the dielectric response of all the media involved, the 
vdW-Ld interaction can have both signs^^. 

One also notes that the contributions of the radial and the 
axial directions are inherently coupled in the Hamaker coef- 
ficients. This means that we cannot simply average the in- 
teractions of the pure axial £ e || and radial e c j_ components, 
particularly if there is a large degree of optical anisotropy in 
a given material. Optical contrast is thus the root source for 
vdW-Ld interactions. Geometry determines how that contrast 
is weighted, while the configuration, or stacking order of the 
materials and thereby their vdW-Ld spectra in the calculation, 
determines the sign. There are additional complexities (e.g. 
retardation^, presence of many walls in the CNTs, surfac- 
tants covering the CNTs^ etc.), but these topics are beyond 
the scope of this paper and have been treated at length else- 
where^. 

6 van der Waals - London dispersion spectra 

Of the three types of e"(co) variations: position, shape and area 
of the peak, SWCNT vdW-Ld spectra are primarily affected 
by shifts in Jt] — JC5 dispersion peaks as well as the details of 
the type of metal peak. In the first part of this paper we defined 
and outlined these differences in the electronic structure and 
dielectric response functions in detail. Here we limit ourselves 
only to the most noteworthy vdW-Ld spectra. 

6.1 Armchair SWCNT vdW - Ld spectra 

VdW - Ld spectra e(z^) of armchair SWCNTs are the most 
straightforward class for several reasons: 1) they have identi- 
cal cutting line or chiral angles due to their similar geometrical 
construction, 2) they are all electronic structure metals, 3) they 
are all null metals, and 4) pragmatically they are by far the eas- 
iest to calculate. Points 1, 2, and 3 are advantageous. We can 
focus solely on the effects of radius before adding additional 
features (e.g. metal peaks, metal/semiconductor differences in 
the Jti — TC5 peak positions, etc). Point 4 becomes an issue for 
chiral SWCNTs, which require a larger number of atoms in the 
band structure calculation and therefore more computation. 

Large diameter armchair SWCNTs have nearly identical a 
peaks, no metal peaks, and Jti — 715 peaks that shift as a func- 
tion of inverse diameter. Thus the resulting vdW-Lds are pri- 
marily affected by the dispersion properties in the 0.1-5 eV 
range. Figure [8] compares e"(co) and for the armchair 
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axial and radial directions of several large diameter SWCNTs. 
e"(co) and the resulting e(;cj) show smooth and systematic 
trends. Note the different variation in the high/low energy 
wings of the dispersion spectrum with respect to radius. While 
the variation of the radial part of the dispersion spectrum is ob- 
viously very different from the variation of the axial response, 
the final effect on the vdW-Ld spectra is very similar. Both 
show a general increase in the magnitude of the low energy 
wing with increasing radius. 

There is a clear difference in the behavior of the radial and 
the axial component of the dispersion spectrum for a given 
SWCNT, Figure [8] In our discussion of the optical proper- 
ties above we outlined the notable difference in the 7Ci posi- 
tion between the electronic structure metals and semiconduc- 
tors. Briefly, the first and subsequent K peaks for an electronic 
structure metal typically occur at a much larger energy than 
in a semiconductor of equivalent radius. The same is true be- 
tween the axial and radial directions of the armchair SWCNTs, 
which behave as electronic structure metals in the axial direc- 
tion and semiconductors in the radial direction. The result is 
that semiconducting SWCNTs have their first % peak much 
closer to eV, resulting in a steeper vdW-Ld spectrum around 
eV. 

As shown above, the regularity in the trends of dispersion 
peaks with respect to radius began to break down at SWCNT 
radii of approximately 0.7 nm (e.g. the zone-folding limit). 
However this breakdown in behavior of dispersion peaks does 
not entail a corresponding change in the vdw-Ld spectra for 
smaller radii tubes. Small changes in the shape of e"(co) are 
not enough to impact the vdW-Lds. Figure [8]compares the re- 
sulting axial and radial direction vdW-Ld spectra for the arm- 
chair SWCNTs ranging from the [24,24,m] all the way down 
to the [3,3,m]. These results show that the trends are system- 
atic down to [5,5,m], which has a radius of 0.34 nm or about 
half the breakdown diameter size for the dispersion peaks. 
Therefore, we can expect that other SWCNT classes will ex- 
hibit similar resiliency of the systematic trends in the behav- 
ior of vdW-Lds for smaller diameter SWCNTs. However, 
SWCNT types that are less symmetric (i.e. chiral semimetals 
and semiconductors) are less likely to have such clean trends; 
the larger variation in their interband transition peaks in the 
dielectric function is caused by a more diverse set of atomic 
structures. 



6.2 Metal SWCNT Ld spectra 

Metal peaks can contribute strongly to vdW-Ld spectra, 
Hamaker coefficients, and consequently to total vdW-Ld inter- 
action energies. According to the classification introduced in 
the first part of this paper, we group the three types of behavior 
in the dielectric function that were observed in the 0.00-0.10 
eV metal regime into null, M0, and Ml classes. Despite this, 



SWCNTs having a continuous DOS through the Fermi level, 
the armchair SWCNTs have no metal features while the chi- 
ral metals exhibit both M0 and Ml metal type peaks. This 
is the metal paradox whereby symmetry forbids certain inter- 
band transitions, giving them a vanishing dipole matrix ele- 
ment despite the available conduction to valence band states. 

The divergence of axial vdW-Ld spectra near t, — depends 
strongly on both the size and location of the interband transi- 
tion peaks. Thus, for a given peak magnitude and area un- 
der the curve, an M0 peak would result in a larger vdW-Ld 
spectral amplitudes since its peak in the dielectric function is 
shifted to eV. However, an Ml peak can still have a larger 
effect than a M0 peak if its oscillator strenght is significantly 
greater. One thus cannot directly correlate one dielectric func- 
tion peak with one type of vdW-Ld spectral feature because 
the integral transform between the dielectric function and the 
vdW-Ld spectrum. Yet there is a tendency for the metal M0 
peaks to contribute more strongly for a given peak height and 
area under the peak. 

The M0 and Ml axial peaks can sometimes vary systemat- 
ically as a function of radius within a given SWCNT type. 
In the first part of this paper, the zigzag type of SWCNTs 
are found to have M0 peaks for the [9,0,m] all the way to 
the [21,0,m], varying inversely with radius. For zigzag met- 
als larger than the [21,0,m], the M0 peaks give rise to a Ml 
peak with a magnitude that eventually became negligible by 
the [30,0,m] SWCNT, Figure [9] 
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Fig. 9 Axial e(i^) of zig-zag metals for energies close to eV. The 
vdW-Ld spectrum e(i^) is obtained from the appropriate e"(co) 
presented in Fig. [8] Frequency in eV. 

In comparison, chiral metals lack the systematic shift in M0 
and Ml peaks as a function of radius; consequently their vdW- 
Lds magnitudes are harder to classify. Additionally, the differ- 
ences in their o and K peak positions as a function of their chi- 
ral angles result in further changes of spectral features across 
the entire energy range. Thus the armchair, zigzag, and chiral 



metals can offer very diverse properties. Figure 10 shows spe 
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Fig. 8 A comparison of radial and axial e"(co) (top) and vdW-Ld e(i^) (bottom) variations in the radial and axial directions for armchair 
SWCNTs, from [15,15,m] to [24,24,m]. Radii and electronic structure (metal, semimetal, semiconductor) of the SWCNTs are listed in Table 
1. Frequency in energy units (eV). 



cific examples of radial and axial vdW-Ld spectra for similar 
diameter SWCNTs within each structural class. 



6.3 Semiconductor SWCNT Ld spectra 

Large-diameter zigzag semiconductor SWCNTs are the last of 
the three classes of SWCNTs in which we can observe the ef- 
fects of increasing radius while fixing the chiral angle and thus 
limiting variability in the o interband transition peaks. Obtain- 
ing optical properties such as the dielectric function for very 
large diameter chiral SWCNTs is still difficult but will likely 
become easier with greater computational resources. Figure 
1 1 compares vdW- Ld spectra for large zigzag semiconductor 
SWCNTs with those of the armchair metal SWCNTs. Be- 
cause these SWCNTs are semiconducting in the axial direc- 
tion, their Tti peak occurs much closer to OeV and results in 
a sharper vdW-Lds wing near eV. Because of the lack of a 
metal peak and a shift of Tti peak out to higher energies, the 
armchair metals have a vdW-Lds slope that is comparatively 
much flatter. However, the behaviors of the radial vdW-Ld 
spectra for these two geometrical types are similar because 



their % dispersion properties are largely unaffected by the rela- 
tionship of the cutting lines with the K-point crossings. Varia- 
tions in the o properties are identifiable but may be ineffective 
as a change in shape versus a change in dispersion spectrum 
peak-width or position. 

In the armchair classification, variations in the dielectric 
function appeared to cancel and result in smoothly varying 
vdW-Lds down to the smallest SWCNTs. However, small di- 
ameter chiral SWCNTs were found to have considerably more 
spectral variation because of the changing electronic structure 
and cutting line angles. The result is a collection of vdW-Lds 
that don't appear to behave in any systematic or definable way. 
Figure 13 compares e"(co) and of armchair and chiral 



SWCNTs of varying electronic structure properties. Any cor- 
relations are inconsistent and unpredictable. The [6,2, s] and 
[6,4,s] have nearly identical axial direction e(/cj), but differ 
greatly from [6,1, s], which has an axial e(/^) similar to that 
of [5,1, s]. While a given pair may have similar axial 
they might still have completely different radial e(/^), and vice 
versa. 
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Fig. 10 The radial (left) and axial (right) vdW-Ld spectra of SWCNTs representing all three structural types of electronic structure metals: 
armchair, zig-zag, and chiral. Examples of radial and axial vdW-Ld spectra for similar diameter SWCNTs within each structural class. For 
details see Table 1. Armchair SWCNT [10,10,m], [1 1,1 l,m], zigzag SWCNT [15,0,m], [18,0,m], and chiral SWCNT [8,5,m], [9,3,m], 
Frequency in eV. 




Energy (eV) Energy (eV) 



Fig. 11 The axial (left) and radial (right) vdW-Ld spectra for the armchair ([23,23,m], [24,24,m]) and electronic structure semiconducting 
zig-zag ([28,0,m], [29,0,m]) SWCNTs. Radii and electronic structures from Table 1. Note that the vdW.Ld spectra of different types of 
SWCNTs can cross. Frequency in eV. 



7 Hamaker Coefficients 

Do these trends in the electronic structure and corresponding 
vdW-Ld spectra persist to the level of the Hamaker coefficient 
and vdW interactions? Figure (a) shows the + Sl^ 
Hamaker coefficients derived in Sec. 5.2 at small separations 
for cylinder-cylinder interactions for all 64 SWCNTs across 
vacuum and water, using Eq. [5] Intricate and systematic vari- 
ations motivate further examination. 

In the limit of large CNT radii, the Hamaker coefficients for 
different classifications of SWCNTs trail different asymptotes, 



Figure 12 The numerical values of the Hamaker coefficients 



are relatively flat and unchanging above the zone-folding limit 
corresponding to a CNT radius of about 0.7 nm. Below this 
limit, their magnitudes begin to shift upward as an inverse 
function of diameter. One observes a widening envelope that 



encompasses all the variations in the Hamaker coefficients and 
becomes larger as the tube radii get smaller. For the smallest 
tubes, the difference in the Hamaker coefficients can be almost 
an order of magnitude between different classification types. 

The large variation in vdW-Ld spectra as a function of (n,m) 
results in a substantial variation in the corresponding Hamaker 
coefficients, a result that could not have been predicted by the 
tight binding approximation, electronic structures classifica- 
tion analysis, force-field modeling, or other approaches other 
than the explicit determination of the electronic structure and 
optical properties and the resulting full spectral Hamaker co- 
efficients. 

Hamaker coefficients and vdW-Ld interactions show a wide 
variation among different CNT types. To simplify we limit 
the discussion here to the particular exemple of two p aralle l 
cylindrical SWCNTs at the small separation limit (se d 34 * 35 !)- 
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Fig. 12 Comparison of Hamaker coefficients, Eqs. [5] vs. SWCNT 
radius a(nm) for a cylinder-cylinder geometry at small separations, 
in vacuum and water. Top: All calculated SWCNTs in vacuum vs. 
in water. Bottom: Comparing the 5 vdW-Ld classifications in vacuo. 
The water dispersion spectrum is canonical, taken from Ref.1251 



A system design example will be included in the remain- 
ing two subsections to demonstrate practical usage (See also 
the Gecko Hamaker project, |http: //sourcef orge . net 7] 



pro jects/geckopro a free and open source project that en- 
ables the calculation of full spectral Hamaker coefficients for 
these and many other systems). 

7.1 Hamaker Coefficients and SWCNT Properties 

Figure 



12 



(b) compare Hamaker coefficients for various 
classes: armchair, chiral metals, chiral semimetals, zigzag 
metals and zigzag semimetals. The Hamaker coefficient be- 
haves systematically across all individual classes but exhibits 
substantial variation between classes. This is not surprising 
given the variety of metal and K dispersion peak behavior 
among the different classes. The larger diameter tubes in Fig- 
ure 12 (b) have a clear ranking in Hamaker coefficient strength 
as a function of spectral class. However, these trends begin to 
cross over below the 0.7 nm zone folding limit, and there is 
also some additional noise/variation from the introduction of 



the chiral classes. Overall, Hamaker coefficient strengths for 
the SWCNTs drop by a factor of 2 when the radius decreases 
from the 0.2 nm to 0.5 nm. For tubes larger than 0.5 nm, the 
decrease in Hamaker coefficient continues through the largest 
tubes examined here. 

We cannot completely attribute the larger range of Hamaker 



coefficient magnitudes in Figure 12 to the effects of increasing 
curvature at the small diameter limit. The large diameter limit 
represents only 3 of the 5 SWCNT categories because of the 
computational difficulty in obtaining very large diameter chi- 
ral metals and semiconductors. It is possible that the variation 
at this limit will be seen once these last two types are included 
at this diameter limit. 

The magnitudes of Hamaker coefficients for larger radius 
SWCNTs in no way prevents the ability to separate these 
SWCNTs. To separate by class, one would design meth- 
ods that target the different asymptotic limit behavior of the 
Hamaker coefficients. To separate within a class (perhaps by 
radius) one would have to create a Hamaker coefficient that 
would divide along the classification, see below. 

7.2 Hamaker Coefficients and the Interaction Medium 



In Figures 12 a) we show the corresponding Hamaker coeffi- 
cients in water and in vacuo. For all but the smallest SWCNTs, 
the Hamaker coefficients differ systematically by a factor of 2 
between these two media. 

7.3 Hamaker Coefficients and SWCNT Size 

For macroscopic systems, material properties are typically in- 
dependent of size. However, for small scales the material 
properties and size can be coupled, as is the case here for 
SWCNTs whose dispersion spectra are connected with their 
radii. The same is true at very large diameters, where vdW-Ld 
spectra and Hamaker coefficients can be again very similar. 

For example, the 1% drop in the Hamaker coefficient be- 
tween the [18,18,m] and [24,24,m] SWCNTs is small rela- 
tive to the 33% increase in the radius a. However, the vdW- 
Ld interaction energy varies as -^/a, Eq. |4] This means the 
[18,18,m] and [24,24,m] may have a Hamaker coefficient vari- 
ation of only 1%, but the total energy varies by 13%. In short, 
Hamaker coefficients can act as useful guides to gauge the in- 
teraction strength but they cannot be the sole source for cal- 
culating the total vdW-Ld interaction energy that contains ex- 
plicit dependence on the radius of the tube. 

Figure 14 shows Hamaker coefficient multiplied by y/a. 
We see two competing effects. As the Hamaker coefficients 
themselves drop as a function of radius and reach an asymp- 
totic behavior around 0.8 nm radius and above, variation of 
the Hamaker coefficient multiplied by ^fa, i.e., the total inter- 
action free energy, for each SWCNT type results in a bowl- 
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between two SWCNTs can thus be varied either by chang- 
ing the materials, affecting the Hamaker coeffcient directly, or 
by changing the radius of the SWCNT. However, the SWCNT 
material properties and their radii are coupled. From Figure 14 
one can deduce that in the small-radius regime changes in the 
interaction free energy at a given surface separation are more 
directly correlated with the changes in the Hamaker coefficient 
itself, while for larger radii there is a small difference in the 
Hamaker coefficients but the changes in the radius can be big. 
We term the first type of behavior "material dominated" and 
the second one "size dominated". 



Fig. 14 Comparison of Hamaker coefficients multiplied by ^fa for 
two rodlike SWCNT cylinders in the limit of small separations in 
water. Hamaker coefficients were calculated using the anisotropic 
rod-rod close separation result with armchair vdW-Ld spectra. The 
interaction free energy, proportional to Hamaker coefficients 
multiplied by ^fa at a fixed spacing, distinguishes "material 
dominated" and "size dominated" behavior for small vs. large radii. 



shaped curve. The interaction free energy at a fixed separation 



8 Conclusions 

We have reviewed e"(ca) optical dispersion properties for tran- 
sition energies up to 30 eV and Hamaker coefficients of the 
vdW - Ld interactions of 64 SWCNTs. The resulting trends 
in the 0-5 eV or 7t regime were found to depend on the cut- 
ting line position and packing, while the 5-30 eV or a regime 
depended largely on SWCNT geometry. We discovered that 
the trends as a function of radius tend to break down and be- 
come more erratic for tubes of radius less than 0.7 nm. The 
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trends, breakdowns, and other interesting results can be cat- 
egorized in a new SWCNT classification based on a combi- 
nation of electronic structure, curvature, and geometry. Ulti- 
mately, the analysis here can lead to a further understanding 
of the chirality-dependent properties of SWCNTs exploitable 
by various experimental means. The use of ab initio methods 
to obtain realistic directionally dependent optical absorption 
spectra for highly complex systems, that are not available by 
other means'^, can accelerate the study of vdW-Ld interac- 
tions in these systems. 

We discovered that the optical properties for SWCNTs need 
a new optical dispersion classification scheme in order to ade- 
quately group and address certain features that would be perti- 
nent to practical properties. We have classified the sub-groups 
of SWCNTs with regard to a, K and metal features. These 
properties determine also the Hamaker coefficients and total 
energies of vdW-Ld interactions in many ways to serve as a 
basis for novel classification schemes. Standard classifications 
based on fewer groupings of the properties (e.g. zone-folding 
and Lambin classifications SJ) or on symmetry considerations 
(e.g. Saito/Dresselhaus classification 9 ) would in our opinion 
neglect key possibilities in experimental design. Exploiting 
important differences in the vdW-Ld interactions of nanoscale 
objects will allow better separation and placement into nan- 
odevices. 

Some SWCNT classes may ultimately be redundant (e.g. 
the chiral zigzag semiconductors) because the differences are 
not significant enough to warrant such fine distinction. Our 
optical dispersion classification is also equally incapable of 
predicting the unique features that occur for the very smallest 
SWCNTs. For example, it cannot explain why the [5,0,m] is 
a metal despite having an (n,m) vector that predicts a semi- 
conductor. When curvature effects become large, the impact 
on e"(co) becomes too strong to allow for a generalized cat- 
egorization. The optical dispersion classification also fails to 
explain why the large diameter zigzag metals lose their metal 
features. 

There are also other material properties that may require 
revisiting this classification system. The effects of structural 
relaxation^, phonons^, bending moments^, functionaliza- 
tion^, the importance of excitonic effects in the optical spec- 
tral calculation 5 1 " 55 > and interlayer coupling^ may in fact 
cause armchair tubes to lose some symmetry and regain a low- 
energy spike. Furthermore a high tran sverse ele ctric field may 
in fact cause the band gaps to open U pl2il23121187j Th^g effects 
are well beyond the scope of this paper, but should be kept in 
mind in order to see if some differentiations are redundant or 
if new sub-classes need to be formed. Even if changes in the 
metal and n peaks do arise, the more stable a peaks should 
still justify the differentiation from at least an optical property 
characterization standpoint. Features of SWCNT dielectric re- 
sponse justify their newly proposed dispersion classification. 



While most of this analysis was specific to SWCNTs, the gen- 
eral concepts described can be applied to other nano-materials 
where vdW-Ld spectral properties can be altered as a function 
of geometry or other controllable parameters. 
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Table 1 Tabulated data of all SWCNTs. The chirality vector (n,m) 
determines the radius, angle, and geometry. The classification 
depends on the relationship between (n,m) and symmetry. The 
number of atoms required for a lattice repeat in the ab initio 
calculations is also included. 
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Fig. 15 Ab initio One-electron theory based calculations to obtain the band structures is used to obtain the frequency dependent dielectric 
response for 64 SWCNTs and based on this to compute their Van der Waals - London dispersion interaction energies. 



22 | Journal Name, 2010, [vol], 1-[22] 



This journal is©The Royal Society of Chemistry [year] 



